Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 



(MN MfeK style file v2.2) 



Halo Concentration and the Dark Matter Power Spectrum 



Kevin M. Huffenberger & Uros Seljak 

Department of Physics, Jadwin Hall 
Princeton University, Princeton, NJ 08544 



January 2003 



o 
o 

(N 

c 

a 



m 

o : 

m . 

o ■ 

6 ■ 
a: 



13 



ABSTRACT 

We explore the connection between halo concentration and the dark matter power 
spectrum using the halo model. We fit halo model parameters to non-linear power 
spectra over a large range of cosmological models. We find that the non-linear evolu- 
tion of the power spectrum generically prefers the concentration at non-linear mass 
scale to decrease with the effective slope of the linear power spectrum, in agreement 
with the direct analysis of the halo structure in different cosmological models. Using 
these analyses, we compute the predictions for non-linear power spectrum beyond the 
current resolution of iV-body simulations. We find that the halo model predictions 
are generically below the analytical non- linear models, suggesting that the latter may 
overestimate the amount of power on small scales. 



1 INTRODUCTION 



In recent years, computational iV-body simulations have 
provided important insights i nto the non-linear form ation of 
structure by the dark matter. INavarro et alJ Jl997f) (NFW) 
observed a 'universal' density profile for dark matter haloes, 
suggesting that the profiles depend only on the mass of the 
halo, with relatively little scatter between them. They also 
noticed that the halo concentration parameter, describing 
the degree of profile steepness (defined more precisely be- 
low) decreases with increasing halo mass, though the exact 
relation depends on the cosmology at h and. Several formulae 
have been suggested for this relatio n iNavarro et al.lll997t 
iBullock et alJl200lt lEke et ablkoOll) which rely on fitting 
some mechanism to haloes culled from iV-body simulations. 
These studies found that the concentration dependence on 
the power spectrum is caused primarily by the slope of the 
linear power spectrum. Both this and the mass dependence 
of concentration can be explained by the epoch of the halo 
formation, since less massive haloes form at higher redshift 
when the universe was denser and this leads to a more con- 
centrated profile. 

Another relatio n which appears to be universal is that of 
the mass function II Jenkins et alJl20011) . The distribution of 
halo masses in A^-body simulations appears to be universal 
in the sense that when mass is related to the fluctuation 
amplitude of the linear power spectrum, the mass function 



detail in recent work 


ISheth & Tormenllf 999l; iJenkins et alj 


200 111 White! 12001. 


2002), defining precisely what is meant 



by the mass and confirming the universality of the mass 
function for a large set of cosmological models. 

While the above approaches have focused on individ- 
ual haloes in the simulations by counting them and explor- 
ing their structure, for many cosmological applications all 
that we need is the dark matter power spectrum. This can 
be computed from A^-body simulations as well, but com- 



putational limitations prevent one from exploring a large 
range of cosmological models or from achieving a high reso- 
lution on small scales. For this reason many fitting formulae 
have been developed w i th increasing accuracy over the years 
jHamilton et alJ 119911 IJain et alJ If 9951 : iPeacock fc Doddsl 
1996; Smith et al. 2002), although these are still significantly 
limited on small scales by the dynamical range of simula- 
tions. 

The close connection between the two descriptions 
has been highlighte d recently with the revival of the halo 
model <SeliaJj200d|Peacock fc Smitrl200d:lMa fc Frvl200(i 
IScoccimarro et alJl200ll) . For a few cosmological models, it 
was shown that with an appropriate choice of concentra- 
tion mass dependence and mass function one can use the 
halo model to accurately predict the non-linear dark mat- 
ter power spectrum. One of the questions we explore in this 
paper is whether this agreement can be extended to a wider 
range of models and whether the trends seen in the analy- 
sis of individual haloes are confirmed also with the power 
spectrum analysis. There are several reasons why this is of 
interest: since the selection of individual haloes is somewhat 
subjective it is possible that those left out may be system- 
atically different. For example, they could be less relaxed, 
more recently formed and so less concentrated. An opposite 
effect is that because the power spectrum is a pair-weighted 
statistic, if there is a scatter in the mass-concentration re- 
lation then pair weighting increases the mean concentration 
relative to the simple particle weighting. 

On the other hand, if the analysis of individual haloes 
is in agreement with the power spectrum analysis over the 
range where both are reliable, then the same analysis can be 
extended to smaller scales which are not resolved by cosmo- 
logical Af-body simulations that compute the power spec- 
trum. The reason for this is that the resolution scale for the 
halo structure can be extended significantly by resimulating 
the representative regions of the haloes with a higher resolu- 
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tion simulation zoomed on that halo. On the other hand, a 
power spectrum calculation requires a large simulation vol- 
ume, so that the largest scales are in the linear regime, which 
then implies that the mass and force resolution cannot be 
as high as in the simulations which focus on single haloes. 

This paper is structured in the following way. We sum- 
marize the halo model and give relevant background infor- 
mation in In JJ] we give our method and results. In 2] we 
compare our results to other models and in JSlwe summarize 
our conclusions. 



2 BACKGROUND 
2.1 The Halo Model 

We begin by reviewing the formalism of the halo model to 
compute the no n- linear dark matter power spectra, following 
the notation in Seliak (2000). The haloes are characterized 
by their mass M, density profile p(r,M), number density 
n(M) and bias b(M). iV-body simulations suggest a family 
of density profiles: 



p(r) = 



(1) 



(r/r s )- a (l +r/r s ) 3 + a ' 
where p s is a characteristic density, r s is the radius where the 



profile has an effective l 


jower 


law index of —2, and —1.5 < 


a < — 1 (Navarro ct al. 


1997: Moore et al.lll998fl. Here we 



use a — —1, since power spectra are not sensitive to the 
inner parts of the halo. 

As is conventional, we re-parametrize p s and r s in equa- 
tion Q in terms of a halo mass and concentration. We 
define the virial radius r v i r of a halo to be the radius of 
a sphere with some characteristic mean density, discussed 
below. Insisting that the mass contained inside r v j r is M 
fixes p s for a given r s . We introduce the concentration pa- 
rameter, the ratio c = r v it/r s . There are several defini- 
tions of virial radius used in the literature. In this paper, 
the virial radius is defined as the radius of a halo-centreed 
sphere which has a mean density 180 times the mean den- 
sity of the universe. This may be denoted rigoo. Other au- 
thors set r v i r to the radius inside which the mean density is 
200 ti mes the critical dens ity (r2oo), a measure independent 
of dNavarro et al.l fl997l . for example). Still others use a 
spherical collapse model to estimate the radius of a vital- 
ized halo (t-a). These lead to different concentrations and 
masses. For f2 = 0.3 and A = 0.7, r2oo < ?"a < nson, so for 
the same halo with Cison = 10 one has C200 ~ 6, ca ~ 8, 
M 2 oo/Mi8on ~ 0.72 and Ma/Misoq ~ 0.87. Assuming the 
NFW halo profile it is straightforward to translate values 
for concentrations and masses between these conventions. 
We choos e cison as it permits the use of the universal mass 
function J Whitel 12002k To avoid cumbersome notation, in 
this paper we refer to Cison as c and Mison as M. 

For the mass dependence of concentration we choose a 
power law parametrization 



c = c (M/M*y 



(2) 



where Co and /3 are free parameters, and M* is the non-linear 
mass scale, which we will define shortly. This is of course not 
the most general parametrization, but as we will show below 
it suffices for the current dynamical range. 



The halo number density is written in terms of the mul- 
tiplicity function via 



n(M) dM = jLf(y) dv, 



(3) 



where p is the mean matter density of the universe and mass 
is written in terms of peak height 



S c (z) 
<t(M) J ' 



(4) 



where S c (z) is the spherical over-density which collapses at 
z (5 C ~ 1.68) and a(M) is the rms fluctuation in the matter 
density smoothed on a scale R = (3M/47rp) 1/3 . v(M,) = 1 
defines the non- linear mass scale. A cosmology with fio = 
0.3, A = 0.7, T = 0.2, and a 8 = 0.9 has M, » 10 13 ft" 1 M e . 

The power spectrum has two terms. The first corre- 
sponds to correlations in density between pairs of points 
where each member of the pair lies in a different halo, and 
so is named the "halo-halo" (hh) term. The second corre- 
sponds to correlations between pairs in the same halo, and is 
known as the one-halo or Poisson (P) term. For convenience 
in calculating convolutions, we work in Fourier space, intro- 
ducing the Fourier transform of the halo profile, normalized 
by the virial mass, 



y(M,k) 



1 

M 



2 , . sin(fcr) 
47rr p{r)- 



kr 



dr. 



(5) 



The mass of the NFW profile is logarithmically diver- 
gent, so in order to evaluate this integral, we must impose 
a cutoff. This does not have to be at the virial radius, since 
we know that haloes are not completely truncated there. In- 
stead, NFW profile typically continues to 2 — 3r v i r . In this 
regime there is already some overlap between the haloes, so 
for the purpose of correlations one can count the same mass 
element in more than one halo. Thus the mass function inte- 
grated over all the haloes may even exceed the mean density 
of the universe (but it can also be below it, since it is not 
required that all the matter should be inside a halo). 

The halo-halo contribution is 



P hh (k) 



Plin(fc) 



dvf(v)dv b{y)y[M{y),k) 



(6) 



where b(v) is the (linear) bias of a halo of mass M ( v), fo r 
which we use the A^-body fit of ISheth fc Tormenl <199Sft . 
Since we want this term to reproduce the linear power spec- 
trum on large scales w e impose this constraint onto the form 
for b(v) |Scliak 200(3). The Poisson contribution is 



duf{u) MM\ y{M (v),k)\ 2 . 



(7) 



If a cutoff beyond the virial radius is used, the mass weight- 
ing should reflect the increase. The treatment of the Poisson 
term on large scales is only approximate. Mass and momen- 
tum conservation require that on very large scales the non- 
linear term should scale as k 4 , rather than as a constant 
implied by equation |7] so the contribution from this term 
on large scales is overestimated for k r~ ; J. This partially 
compensates the increase in power from the matter outside 
the virial radius and for this reason we chose to use r v i r as 
the radial cutoff for the halo. The halo-halo term is only ap- 
proximate, since we do not include the exclusion of haloes, 
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Figure 1. Using the halo model, several power spectra were gen- 
erated from a cosmology with f2o = 0.3, A = 0.7, T = 0.2, and 
erg = 0.9 at z = 0. In the upper panel, cq varies and f3 = —0.2. In 
the lower panel, /3 varies and co = 10. 

which would suppress the term. Because of these approxi- 
mations one would not expect the halo model to be perfect, 
especially in the transition between the linear and non-linear 
regime. 

The total power spectrum is the sum of the two contri- 
butions, P(k) = P hh (k)+P p (k). On small scales, the Pois- 
son term dominates. On larger scales, the halo-halo term 
dominates and reduces to the linear power spectrum. From 
here forward, we express the power spectrum in dimension- 
less form, A 2 (fc) = Ank 3 P(k). 

The effect on the power spectrum of varying Co and B in 
a typical ACDM cosmology is shown in Figure The con- 
centration controls how tightly matter is correlated within a 
single halo. Therefore, higher concentration means a larger 
one-halo term. Increasing Co increases the concentration at 
the non-linear mass, affecting the amplitude of the power 
spectrum where it is dominated by one-halo term. 

Varying 3 keeping Co constant produces a tilt in the 
non-linear power spectrum around the fiducial value, which 
is set by the scale where the non-linear mass dominates the 
power spectrum contribution (k ~ 30-40 h Mpc -1 ). Haloes 
less massive than the non-linear mass dominate at higher 
k. Steeper (more negative) means that haloes less mas- 
sive than the non-linear mass will have their concentrations 



enhanced, leading to an enhanced one-halo term, and more 
power at high k. Haloes more massive than the non-linear 
mass will have their concentration reduced. These dominate 
at intermediate k in the non-linear regime, so power there is 
reduced. The changes in shown in the figure may not have 
much effect on the power spectrum at k < 100 h Mpc -1 , 
which is the range where Af-body simulations are reliable. 
On the other hand, they have a substantial effect on the 
concentration of typical-sized haloes. For the same concen- 
tration at the non-linear mass (~ 10 13 /i -1 Mq), a halo at 
10 12 /i -1 Mq has a concentration 50% lower with = —0.1 
than with = —0.4. These masses do not dominate the 
power spectrum until k > 100 h Mpc -1 , so this change does 
not make as much of a difference on the power spectrum at 
scales larger than this. This discussion suggests that while 
power spectrum analysis cannot provide strong constraints 
on the mass dependence of concentration, one can use con- 
centration mass relation to predict the power spectrum at 
small scales which are not resolved by simulations. 



3 METHOD AND RESULTS 

We consider two sets of cases which have been extensively 
simulated, self-similar initial conditions and more realistic 
cold dark matter initial conditions. In both cas es we use the 
fitting formulae for non-linear A 2 (fc) given bv lSmith et alJ 
J2002li . 

3.1 Self-Similar Case 

The simplest case to consider is fl = 1 Einstein-de Sitter 
universe with a power law linear power spectrum. In this 
case, has an analytic form, provided we assume that once 
a halo collapses, the scale radius r s is fixed in proper coordi- 
nates. This is suggested by the simple model which assumes 
the haloes remain unchanged once formed, which seems to 
hold in numerical simulations feullock et al.ll200ll) . We fol- 
low the evolution of a single halo, as it traces out a portion 
of the c(M) relation. 

In an fi = 1 universe, density perturbations grow as the 
scale factor a = 1/(1+ z). For a power law linear spectrum, 
this means Pn n (k) oc a 2 k". Top-hat smoothing at the scale 

_ Tl + 3 

corresponding to M yields a(M) oc aM b . For an Q m = 1 
universe, 8 C ~ 1.68 is constant in time, so a(M.) is also 

constant. This means M„ 6 oc a. If r s is fixed in time, 
c oc rvi r oc a follows from the definitions of c and rvir and 
from p oc a~ 3 . Since c oc MZ^ then (3 ~ —^^g 1 ^ if we assume 
M is constant. In reality M also increases with a, which 
decreases /3 somewhat, although not by more than 20%. 

We calculated halo model power spectra with n = —2.0 
and n = —1.5, for which /3 = —0.16 and —0.25 as calculated 
above, with several values of Co. These spectra are shown 
in Figure |21 The agreement is quite good, given all the lim- 
itations of the halo model. At higher k in the n = —2.0 
case, it appears that the slope of the power spectrum dis- 
agrees independent of Co. For k/k, > 2 this model is better 
fit by a power spectrum with 8 = —.037, cp = 2 .6. How- 
ever, according to figure 12 of ISmith et"al] <|2002T) . simula- 
tion data do not exist for k/k* greater than a few tens in 
the n = —2 case, so the discrepancy is not really significant. 
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Figure 2. Non-linear power spectra in scale free models, divided 
by the linear power. The power law linear spectra are Aj? = 
4irk 3 (k/k t ) n , where n = —2 and n = —1.5. In each case the halo 
model spectra for several different c(M) are shown in solid lines, 
and the fitting formula of Smith et al. (2002) is in dashed lines. 
For the n = —2 case we show three models with /3 = . — 1667, 
which is the analytically predicted value from the text. These 
have Co = 3, 4, 5 in increasing amplitude, and fit poorly. We also 
show the best-fitting fi = —.037, cq = 2.6 with a thicker line. For 
n = —1.5, we show the predicted fi = —.25 with co = 6,8, and 
the best-fitting = -.245, c = 7.1. 



In the n — —1.5 case, which is tested against simulations to 
k/k* = 100, the predicted fi — —0.25 is close to the best fit- 
ted value. This agreement therefore confirms the assumption 
that haloes have a fixed scale radius once they formed. 



3.2 Cold Dark Matter Models 

Next we consider several flat CDM models (we do not con- 
sider open or closed models here, since they are observation- 
ally di sfavored). To comp are the halo model power spectrum 
to the lSmith et alJ <|2002T) power spectrum, we use a simple 



X statistic: 



JV 
i=l 



(8) 



We distributed the N = 40 sample ki evenly in log k, consid- 
ering 0.01 < k < 40 h Mpc" 1 . This range takes us from the 
linear regime to some of the highest k that the fi tting formula 
has be en tested with, according to figure 16 of [Smitfi^^ilJ 



2002). For the errors m, we take 30% of the ISmith et alJ 
2002) power. This is somewhat arbitrary, but roughly the 



size of the combined error in the mildly non- linear regime, 
and is probably conservatively large in the fully non-linear 
regime where the Poisson term dominates. We are ignoring 
the correlations between the points, so the actual value of 
X 2 is just a qualitative measure of the goodness of fit. 

For a variety of cosmologies, we minimized x 2 over the 
two-dimensional parameter space of Co and /3 from equa- 
ti qn Et. We employed a Powell minimization as described 
in lPress et alJ |l992) to give a best-fitting Co and fi. All runs 
were given the same initial point: Co = 10 and (3 = —0.2, and 
terminated when the the minimum changed by less than 
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Figure 3. The values of co and (3 at several values of the mass 
density parameter Qo for the fitted halo model power spectrum 
and the models compared in Jl] The vertical line denotes the 
fiducial model with Co = 0.3. The other cosmological parameters 

El 



are given in 



0.1%. We place an error bar on each minimization by calcu- 
lating x 2 on a grid about the minimum point. We compute 
an error contour where exp(— x 2 /2) falls to half-maximum. 
The error bars in Co and f3 show the maximum extent of 
this contour in each direction. With this definition, these er- 
ror bars cannot properly be used to set confidence limits or 
rule out models, but are included simply to gauge in which 
cosmologies parameters are better or worse constrained and 
to explore parameter degeneracies. Indeed, because of the 
range of k we consider, Co and (5 are somewhat degenerate: 
making j3 more negative has a similar effect to lowering Co 
(compare Figure 0. 

We varied one cosmological parameter at a time, choos- 
ing the other parameters from a fiducial model with Qo — 
0.3, Q + A = 1, T = 0.2, a 8 = 0.9, and z = 0. We 
examined the ranges 0.1 < Qo < 1-0, 0.1 < T < 0.4, 
0.6 < as < 1.5, and < z < 1. We plot the variation of 
Co and (3 with Qo with solid lines in Figure [3] The dashed 
lines are discussed in &\ The relation for Co is well fit by: 



c (Qo) = ll(n /0.3)" 



(9) 



while (3 is consistent with being constant around zero or 
slightly negative. Note that the Q dependence depends on 



© 0000 RAS, MNRAS 000, 000-000 



Halo Concentration and the Dark Matter Power Spectrum 5 




-2 -1.5 -2 -1.5 -2 -1.5 



n eS (M t 



Halo model fit I — • — I 

Eke et al. (2001) 

Bullock et al. (2001) 



Figure 4. The values of cq and /3 as a function of the effec- 
tive spectral index at the non-linear scale n e g for the fitted 
halo model power spectrum and for the models compared in 
d The line type indicates which model. We vary one param- 
eter at a time. The fiducial model, given in 43.21 is marked 
the vertical line at n e f£ ~ —1.7. Each parameter increases 
as n e ff increases except for z, which decreases. In their re- 
spective panels, the points (left to right) correspond to: T = 
0.1,0.2,0.3,0.4, a 8 = 0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4, and 
z = 1.0,0.8,0.6,0.4,0.2,0.0. 



the definition of the virial radius and with a different defi- 
nition there would be a different Q dependence. 

Although Co and /3 vary with all of the other cosmologi- 
cal parameters, it is more instructive to present its variation 
with the effective power law index of the linear power spec- 
trum at the non-linear scale, 



iirff 



d(logP Un (fc)) 



d(logfc) 



(10) 



where Af in (fe*) = 1 defines the non- linear scale. The reason 
is that the slope of the linear power spectrum determines the 
epoch of formation of small haloes that merge into the larger 
halo. For lower effective slope n e s these haloes formed later, 
when the density of the universe was smaller, and the final 
halo concentration of the larger halo is also lower. This is 
depicted with solid lines in Figure|l] Our fiducial cosmology 
has n e tf ~ —1.7. Varying £Iq has no effect on n e ff- When 
varying z, a portion of the variation of Co and /3 should be 
due to the change in tl(z). The variation of Co with n e ft is 
consistent with 



Co = ll(fi(*)/0.3)~ 



^V7 



(11) 



while (3 has large errors, but is typically slightly negative. 



4 COMPARISON TO HALO ANALYSIS 



We n ow c ompare our re s ults t o the models of lEke et al] 
feOOll) and lBullock et alJ feOOll) . These are designed to re- 
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Figure 5. A comparison of the variation of the concentration pa- 
rameter with mass for our fiducial ACDM cosmology (jE2|. The 
masses and concentrations here have been translated to a common 
virial radius convention, discussed in the text. The shaded region 
covers power law models within our error contour, discussed in 
the text, and should not be considered a limit on c(M) relations 
of arbitrary shape. 



produce the concentrations of haloes found in JV-body sim- 
ulations. Throughout we translate to our virial convention. 

We begin by considering concentration as a function of 
halo mass for the fiducial ACDM cosmology defined in >I3.2I 
This is plotted in Figure It is clear that our results are 
in agreement with those from direct halo structure analy- 
sis. However, most of the information in our model comes 
from the haloes with masses between 10 13 and 10 14 h' 1 Mq, 
because haloes smaller than this dominate at scales smaller 
than the range we considered for our fit. This means we have 
a relatively short lever arm with which to determine the 
mass concentration relation and consequently we are unable 
to strongly constrain the slope /3. In direct analysis of haloes, 
concentration is a decreasing function of mass, so f3 < 0. Our 
best fitted values for (5 are also typically negative, although 
the errors are large and positive values also cannot be ex- 
cluded. The strongest confirmation of this prediction comes 
from the analysis of scale- free models in previous section, 
where (3 is negative for n = —1.5. It is clear that measuring 
the mass-concentration directly is a better way to determine 
the mass dependence of concentration rather than using the 
power spectrum, particularly at lower masses. 

Since determining mass dependence of concentration 
from the power spectrum does not appear promising let 
us ask the opposite question: how well can we predict the 
power spectrum using the mass concentration relations mea- 
sured directly from the analysis of individual haloes in 
the simulations? Figure [(J shows the power spectra calcu- 
lated with power law approximations to the c(M ) func- 
tions shown in Figure as well as the lSmith et al] J2002T) 
power spectrum. While there is good agrement between 
them for k < 40/iMpc" 1 , the halo models always predicts 
less power at k > 40/iMpc~ 1 compared to the fits to the 
power spectrum. However, the direct fits are not reliable 
in this regime, since the A-body simulations become unre- 
liable for k > 40/iMpc -1 , which is why the fits presented 
above only use information from k < 40/iMpc -1 . 
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Figure 6. A comparison of power spectra generated with different 
c(M) with analytic fitting formulae. Note that we only fit up to 
k = 40 h Mpc" 1 . 



The conclusion from the above comparison is that on 
small scales there is disagreement between analytic formu- 
lae for non-linear power spectrum and the halo model. Since 
the analytic models have not been calibrated with simula- 
tions in this regime it is possible that they are unreliable 
and so the halo model should be us ed instead (indeed, the 
analytic models of lSmith et alJl2002l do not even predict the 
power spectrum for k > 100/iMpc -1 ). Nevertheless, since 
non-linear predictions at high k have been used in the liter- 
ature, espe cially for pred icting the lensing effects on small 
scales fe.g. lMetcalflll999f) . it is important to recognize their 
possible limitations when using them outside the regime of 
applicability. 

It is unlikely that the discrepancy can be resolved by 
modifying the mass function. The existing simulations re- 
solve the mass function down to 10 11 Mq, which is the mass 
that dominates the power spectrum at k ~ 1000/iMpc -1 . 
Substructure could boost the amount of power, since in the 
mass function one counts only isolated haloes, not those 
that are within larger haloes. Subhaloes within haloes con- 
tribute to the clustering in the same way as isolated haloes 
on scales that are comparable to the scale radius. However, 
the abundance of subhaloes is small compared to the abun- 
dance of isolated haloes even at small halo masses, so it is 
unlik ely that this is a significant correction (Bullock et al. 
l200rf) . Similarly, while the scatter in the mass concentration 
relation boosts the amount of power, we find the effect is 
relatively small for the scatter sug gested from simulations 
(Bull ock et al.l i2001; Coorav & Hu 2001) and again cannot 
resolve the discrepancy. 

For Co versus fio and n e ff our fit tolSmith et al] (120021') 
produces very simila r trends to that in lEke^t^dTl200lll and 
iBullock et all (UpOl). This comparison should be reliable, 
since the haloes that dominate are around 10 13 - 10 14 M Q , 
which are abundant and well resolved in the simulations. For 
f3 we also find broad agreement, although the errors from the 
power spectrum method are very large. 



5 CONCLUSIONS 

We have used the halo model and the non-linear evolution of 
d ark matter power sp ectrum (as given by the fitting formula 
of ISmith et"aill2002li to probe the dependence of halo con- 
centration on cosmological parameters. The most important 
parameters for the concentration are the matter density and 
the effective slope of the linear power spectrum at the non- 
linear scale. We present analytical expressions which give 
concetration at the non-linear scale as a function of these 
two parameters. 

We exa mined our result a gain st the models for halo con- 
centration of lEke et all <l200l]) and lBullock et al] (1200 ill . We 
found broad agreement with the trends presented there in 
that we also find that spectra with lower effective slope n e s 
have lower concentration at the non-linear mass. Although 
there is some difference i n the mass -concentration depen - 
dence between lEke et al] <200ll) and IBullock et al] ]200ll) . 
we find no compelling reason to prefer one model to the other 
based on the power spectrum. This is because the dynamic 
range covered by the power spectrum analysis is too narrow 
to obtain a reliable mass dependence of the concentration 
over a wide range of masses. We have found that using these 
models to predict the non-linear power spectrum results in 
a signific antly lower power spectr um on small scales com - 
pared to lPeacock fc Doddsl dl996ft and lSmith et al] fcOOfl . 
This suggests that care must be exercised when using these 
models in the range where they have not been calibrated 
and may indicate that the amount of dark matter power on 
small scales is smaller than predictions from these models. 

For mildly non-linear k, eve n the best fit t ed ha lo mod- 
els show discrepancies from the ISmith et al] (120021) power 
spectrum at the 20-30% level, both for the self-similar and 
C DM models. Some of this discrepancy arises from the fits 
in lSmith et al] ll2002li . which are only accurate at 10% level. 
However, direct co mparison of the halo model to the N- 
body simulations of IJenkins et al] il998T) also reveals dis- 
crepancies at the similar level. This is not surprising given 
the approximate nature of the halo model and argues that it 
cannot be a full replacement for A^-body simulations in the 
era of high precision cosmology. However, using the param- 
eters derived in this paper, the halo model can be used for 
qualitative predictions of the dark matter power spectrum 
over much wider range of models and scales than previously 
available. 
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